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The p\n fractal decomposition: 

Nontrivial partitions of conserved physical quantities 


Vladimir Garcfa-Morales^ 

^Departament de Termodindmica, Universitat de Valencia, E-46100 Burjassot, 

A mathematical method for constructing fractal curves and surfaces, termed the p\n fractal de¬ 
composition, is presented. It allows any function to be split into a finite set of fractal discontinuous 
functions whose sum is equal everywhere to the original function. Thus, the method is specially 
suited for constructing families of fractal objects arising from a conserved physical quantity, the 
decomposition yielding an exact partition of the quantity in question. Most prominent classes of 
examples are provided by Hamiltonians and partition functions of statistical ensembles: By using 
this method, any such function can be decomposed in the ordinary sum of a specified number of 
terms (generally fractal functions), the decomposition being both exact and valid everywhere on 
the domain of the function. 
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I. INTRODUCTION 

A fractal is a geometrical structure with non-integer di¬ 
mension and the property of self-similarity m- Fractals 
abound in present day applications and are being used 
to describe a variety of phenomena in nature, ranging 
from nanoclusters [5] and semiconductors [6] to galaxies 
[7]. Fractals are important in the statistical mechanics 
of complex systems [HHID] and quantum mechanics nn- 
[H]. Mathematically, they can arise from a variety of pro¬ 
cesses: iterated function systems laia, strange attractors 
of chaotic dynamical systems m, critical phenomena 
[HI, cellular automata HHHH], substitution systems |22j , 
bitwise arithmetic [3S], and arithmetic sequences [M] . 
Recently, a construction has been presented in terms 
of the hierarchical structure of a coin division problem 

Hang. 

In this article we present a new construction of a 
wide variety of fractal structures motivated by a gen¬ 
eral mathematical problem of interest to Physics that is 
easy to state and for which we provide a solution: To find 
whether if given an arbitrary Hamiltonian H is it possible 
to split it into the ordinary sum of a given finite number 
A G N o/ (generally discontinuous) functions, bringing 
the Hamiltonian to the following form 

A-l 

n—0 

We call this problem the partitioning of the Hamiltonian. 
With this problem we exclude, of course, the trivial par¬ 
titioning that results from dividing the Hamiltonian H 
by A or by forming any linear combination of rational 
functions of H. We ask whether is it always possible to 
find A nontrivial partitions of H regardless of the form 
of H (of course, the partitions H^ should themselves de¬ 
pend on i? in a nontrivial way). Since we strive to give 
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a general, constructive answer, we must renounce to de¬ 
mand smoothness (continuity and/or differentiability) of 
the parts i7„. 

The answer to this question is, surprisingly, affirmative 
and we provide an explicit, general method, to construct 
each part given H and A. The partition is valid and 
exact everywhere on the energy surface and, by construc¬ 
tion, the HnS generally possess fractal features. Our ap¬ 
proach makes use of a digit function that we have recently 
introduced in a new formulation of classical and quantum 
mechanics m as well as elementary facts of abstract al¬ 
gebra (finite group theory) [^. In the following section 
we first introduce the digit function elucidating some of 
its mathematical properties that we shall then use. We 
explicitly prove that the digit function induces a cyclic 
group structure on a finite set of integers. We then pro¬ 
ceed to establish our main result, which is contained in 
what we term the pXn fractal decomposition: We prove 
that such a decomposition satisfies Eq. 0 thus solving 
the above problem affirmatively. We then extend these 
results to complex-valued quantities and prove a theo¬ 
rem that generalizes the method, significantly extending 
its domain of applicability. We finish with a short dis¬ 
cussion of the mathematical results, sketching their con¬ 
nection to statistical mechanics and showing how they 
constitute a number-theoretic formulation of the latter. 


II. THE DIGIT FUNCTION AND ITS BASIC 
PROPERTIES 

Let p > 1 be a natural number. Every real number x 
can be expanded in powers of p as 

X = qnP^ + aN-iP^~^ + ■ ■ . + ao + a-ip~^ +a-2P~'^ + ■ ■ ■ 

( 2 ) 

where the a^s are the so-called digits of the expansion 
and are either all positive or all negative, with absolute 
values being integers in the interval [0,p — I]. Eq. ([^ is 
the representation of the real number x in base p. We 
shall call p the radix to avoid confusion with other uses 



2 


of the word ’base’ in physics. Eq. ([^ can thus be written 
as 


iiogp klJ [logp klJ 

a; = sign(a;) ^ |afc|/= sign(x) ^ dp{k,\x\)p’' 

k— — oci k——(yci 

( 3 ) 

where the upper bound in the sum gives the total number 
of integer digits N as N = [log^ |x|J +1. Above, the digit 
function dp(fc, |x|) = |afc| of the real number x has been 
defined through the coefficients \ak\- The digit function 
is explicitly given by [27] 


dp{k,x) 


X 


X 

pk ^ 

-P 

pk-\-l 


( 4 ) 


where [.. .J denote the closest lower integer (floor func¬ 
tion) |29l [30]. To see this note that, from Eq. ([^ (we 
consider x nonnegative for simplicity) 


— qnP +a,N-iP + ■ ■ ■ + dk+iP + ctk (5) 


From Eq. ([^ we have, as well, 


= qnP + ttN-iP ^ -b ... -b ttk+ip (6) 


nfc + 1 


with the natural numbers to be presented in a concise 
way. If we concentrate on the real numbers in the unit 
interval [0,1) one can show that, already these, cannot 
be counted with the natural numbers. For let us assume 
an infinite list of such numbers ri, r 2 ,..., ... where 

n G N. Then, the digit expansion in radix p, Eq. (|^ of 
each such number is 

— 1 OO 

= X! P^dp(fc,r„) = ^p“''dp(-A:,r„) (9) 

k— — oo k—1 

Clearly, the number x G [0,1) given by 

OO 

X = ^P~^dp{0,dp{-k,rk) + 1) (10) 

k^l 

is not in the above list since at least one digit of this 
number is different to the ones of each number r^, n = 
1,2... in the list. The existence of this number x proves 
that no bijection between the naturals and the reals in the 
interval [0,1) is possible and, hence, that the cardinal of 
K is larger than that of N. This is the celebrated Cantor’s 
diagonal argument and it finishes the proof. 

Let us now briefly discuss some elementary properties 
of the digit function, that we shall need below: we can 
shift X to X -b np^'^^ (with n being an integer number) 
and the function remains invariant 


whence, by subtracting both equations, we obtain Eq. 

_ 

The digit function dp(fc,x) is plotted in Fig. m for 
p = 5 and k = 0 (continuous curve) and k = —1 (dashed 
curve). It is a staircase of p levels taking discrete integer 
values between 0 and p — 1. Each time that x is divisible 
by the ascent of the staircase is broken and the level 
is set again to zero and a new staircase begins. We note 
that, in Eq. © we have, asymptotically 

lim X = lim 

p—>-oo p—>-oo 

^ sign(x)dp(Llogp|x|J,|x|)pLiogpklJ (7) 

or, if the limit is strictly taken 

lim X = sign(x)doo (0,|x|) (8) 

p—>-oo 

We note that in m X was taken nonnegative throughout 
(and hence |x| = x and sign(x) = I there) and that the 
digit function as defined dp(fc, |x|) in [57| corresponds to 
dp(fc— 1, |x|) here. We make here this notational adjust¬ 
ment to stay more in tune with standard mathematical 
textbooks |3Ij where the radix p expansion is discussed. 
Those textbooks remain at the level of introducing the 
coefficients Ufc’s as in Eq. (|^ but do not discuss the as¬ 
sociated digit function Eq. Q, whose properties we are 
about to exploit here. 

We first note that the digit function allows Cantor’s 
proof of the fact that the real numbers eannot be counted 


Liogp klJ 

sign(x) ^ dp(A:, IxDp'” 


dp(fc, X -b np ■'■ ) = 


X -b np^+^ 


ph. 


X + np^^^ 


r>fc + l 


pK 


P 

np — p 






np 


P 


.P‘ 


|A;+1 


= dp(fc,x) (II) 


i.e. the digit function is p^“'"^-periodic. We have also a 
scaling property that concerns multiplication of x by any 
integer power of p 


dp{k,p^x) = 


p""x 


p'^x 

rr)k-\-\. 


nk—r 


■P 

J?"-"' 

= dp{k — m, x) 


X 


rxk — 7n-\-l 


Since dp{k,x) is an integer 0 < dp(fc,x) < p — I we 


also trivially have 

dp(0, dp(fc, x)) = dp{k, x) - p 


dp{k,x) 


= dp{k,x) 


( 12 ) 

because [dp(fc,x)/pj = 0. Let m G h. By Euclidean 
division, we can always find integers a and b such that 
m = ap + b and b — dp(0, m) and, thus, we also have 

dp(0, n-bdp(0, to)) = dp(0, n-bap-bdp(0, to)) = dp(0,n-bTO) 

(13) 
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where Eq. (Ill has also been used. Another property, 
that follows directly from the definition, is the following 


n = 1 and operate repeatedly with it. We have, from 
Eqs. ( [II| ) and ( [I^ 



= +ndp(fc,^) (14) 

Because of the above properties, it is clear that if 
n G Z is any integer dp(k,n) is an integer in the in¬ 
terval [0,p — 1]. The digit function for k = 0, dp(0,n), 
maps the set of the integers Z to the set of residue classes 
modulo p denoted as Z/pZ: It extracts the remainder of 
n under division by p. Thus the digit function is an ho¬ 
momorphism between both sets. For example, if p = 2, 
d 2 ( 0 , n) = 0 if n is even and dp(0, u) = 1 if n is odd. 

It is a trivial exercise to show that dp(0,n) provides 
indeed a group homomorphism of the integers Z under 
ordinary addition to the finite cyclic group Z/pZ of the 
integers in the interval [0,p — 1] under addition modulo 
p. Let m and n be such integers. The group axioms 
are all satisfied. We have 1) Closure: The digit function 
dp(0, n) yields the remainder of m -|- n under division by 
p which is, of course, an integer G [0,p — 1] as well; 2) 
Associative property: we have 

dp (0, m -I- dp (0, n -I- k)) 

= dp (0, m -I- n -I- fc) = dp (0, dp (0, to -I- n) -|- fc) 

where Eq. (El has been used twice; 3) Neutral element: 
The neutral element is 0, as in the ordinary sum; 4) In¬ 
verse element: The inverse element of to is p — to, since 
dp (0, TO -I- p — to) = 0. 

Thus, the group axioms are satisfied. Furthermore the 
group is abelian since dp (0, to -I- n) = dp (0, n -I- to), and 
it is cyclic because a single element generates the whole 
group. To show this latter statement let us just consider 


dp (0,1 +I) = 2, ..., dp(0,l+p-2) =p-1, 
dp(0,1+p-1) =0, dp(0, l+p) = l 

Thus, the orbit of the element 1 has order p equal to the 
one of the group, 1 can be made the only generator of 
the group, and the group is cyclic. 

The result of operating on to and n with the digit 
function dp(0, to + n) is given in the following table 


dp(0, TO + n) 

0 

1 

2 

... p-2 

p- 1 

0 

0 

1 

2 

... p-2 

p- 1 

1 

1 

2 

3 

... p- 1 

0 

p-2 

p-2 

p- 1 

0 

... p — 4 

p-3 

p- 1 

p- 1 

0 

1 

... p-3 

p-2 


where m is listed in the rows and n in the columns. The 
table gives the output of the operation dp(0, m-\-n) and 
coincides with the Cayley table of the finite cyclic group 
Z/pZ of order p [55]. Specific examples are given by 


d 2 ( 0 ,TO + n) 

0 1 

0 


0 1 

1 


1 0 

d 3 ( 0 ,TO + n) 

0 

1 2 

0 

0 1 2 

1 

1 2 0 

2 

2 0 1 


The above tables have the latin square property char¬ 
acteristic of finite groups and other magmas (as loops 
and quasigroups [32]). Because of the group being cyclic 
each row of outcomes is a cyclic permutation of the pre¬ 
vious row. The abelian character of the group is easily 
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seen through the symmetry of the tables regarding their 
main diagonal. Since all integers S [0,p — 1] appear in 
each row we have the following remarkable fact 

^dp(0,TO + n) = (15) 

n—0 j—0 


This definition means that out of the function f{x) we 
construct the function FpXnf{x) that we shall call a frac¬ 
tal object. There is a total of A such fractal objects, each 
one indexed by n S [0, A — 1]. 

The meaning of this definition shall now be substan¬ 
tiated by proving the following fact, which is the main 
result of this article 


since by fixing m S [0,p — 1] the sum in Eq. (15) runs 
over all possible values n e [0,p — 1] an so do the out¬ 
comes. Note the insensitivity to m of the sum in Eq. 
([T^: The latin square is also a magic square and all 


A-l 


rows and columns add to the same quantity as 

given by Eq. (15). We exploit this property in the next 


section in establishing our main result. 


III. THE pAn FRACTAL DECOMPOSITION 


fix) = FpAn/(a;) 


(18) 


n—0 


i.e. the sum of all A fractal objects yields the function 
from which they constitute a decomposition. Note that 
the function f{x) being decomposed can have any arbi¬ 
trary number of variables f(x, y,. ■.) and any dependence 
on them and the decomposition is equally valid. The re¬ 
sult holds in general. It is straightforward to prove it. 
Note that, because of Eq. (15) we have 


The digit function is a natural building block for fractal 
curves as we shall now explain. In general, as the classical 
example of the Cantor set shows, if we systematically 
exclude from the real line certain intervals, by repeating 
this process at all scales we end up with a fractal set. 
We now describe how a closely related method can be 
mathematically implemented in an efficient manner with 
help of the digit function. 

Let us, for example, set the digits for which |fc| is odd 
in Eq. (|^ to zero. By means of the digit function this 
can be easily accomplished as 

Liogp klJ 

yf{x) = sign(a;) ^ p'=dp(fc, |a;|)d2(0, |fc|-bl) (16) 

k——oo 

this function leads to the fractal curve shown in Fig. 

for p = 3. Note that, when k is odd, the terms 
^dp(/c, |x|)d2(0, |fc| -I- 1) do not contribute to the sum 
because d2(0, fc -I- 1) = 0. In this way, we are systemati¬ 
cally excluding intervals in the real line at all scales. The 
result is a self-similar discontinuous structure (had we all 
the digits we would obtain the line y = x). If we now con¬ 
sider the fractal curve where we set the digits with even k 
to zero yf{x) = sign(a::) p’^^pik, |a:|)d 2 ( 0 , \k\)it 

is clear that the function y{x) = a: is equal to y = yf -\-yf 
since d 2 ( 0 , |fc|) -|- d 2 ( 0 , |A:| -|- 1) = 1. 

We have considered digits in even or odd positions. 
But, having in mind the results in the previous section, 
we can now generalize this as follows. Let f{x), /:!&—>■ 
M be any real-valued function of x (it does not need to be 
neither continuous nor measurable). We define the pXn 
fractal decomposition of f{x) and denote it FpXnf{x) as 

FpXnfix) = (17) 

. . llogpl/(^)IJ 

^ E /dp(A,|/(x)|)d,(0 ,|fc|+n) 

^ ^ k= — oo 


EdA(0,|fc|+n) = ^^^ (19) 

n—0 

Therefore 

^2sign/(x) ^ 

/ . ^P><nf{x) — y ^ ^ 

n—0 n—0 ^ 

Liogp |/(x)|J 

X E /dp(fc, |/(a;)|)dA(0, |fc|-bn) 

k ——oc) 

Liogp |/(a:)|J 

= sign/(a;) E /dp(fc, |/(a:)|) x 

k— — (x> 

Liogp |/(a:)|J 

= sign/(a;) E /dp(fc, |/(a;)|) =/(x) 

k— — oo 

where the radix expansion of f{x), Eq. (§ (i.e. replacing 
X by f(x) in that equation) has been used. 

In Fig. I^two simple examples of this result are shown. 
The functions f{x) = x (left) and f{x) = x^ (right) are 
plotted together with their fractal objects F 3 ^ 3 _„/(a:) (i.e. 
p = 3, A = 3) calculated from Eq. 0 for the values 
of n indicated in the figure. One can see at a glance 
that the sum of the values of the fractal objects is equal 
to f{x) for every value of x. The fractal objects cross 
themselves in nontrivial ways so that their sum is the 
continuous and differentiable function f{x) in the hgures. 
Note that, although the fractal objects are discontinuous, 
they are defined at every point. Since the limits from the 
right and from the left at the discontinuities are different, 
this requires that the jump of the fractal objects at such 
discontinuities agree so as to produce a continuous curve 
when added. This is warranted by the construction. 
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FIG. 2: The fractal curve yf{x) obtained from Eq. (161 for p = 3. 



FIG. 3: (Golor online) The continuous curves f{x) = x (left) and f{x) = x^ (right) and their associated fractal objects 
F 3 , 3 ,nf{x) (i.e. p = 3, A = 3) for the values of n indicated in the figure. 


We have dealt with finite cyclic groups above. In gen¬ 
eral, let g\(rn,n) be a surjective application sending all 
integers m, n to the integers G [0, A — 1]. If the table of 
g\{m,n) (of dimension A^) has the latin square property 
all integers 0,1,... ,p — 1 appear in every column and 
every row. Then, Eq. (15) holds again and we have 


and Eq. 0 can thus be readily generalized as 


pXn 


fix) = 


( 22 ) 


2 sign/(a;) 
A(A-l) 


[logp l/(a:)IJ 


P''dp(fc, |/(x)|)5A(|fc|,n) 


k— — (x> 


A-1 

n—0 


m, n) 


A-l 

= E^' = 

3=0 


A(A-l) 


and the main result Eq. (181 trivially holds since the 


latin square property is all what is needed to prove that 
result. Thus, we can use all finite groups at our disposal 
(21) to generate fractal functions in a nontrivial but system¬ 
atic way, provided that we have a construction for them 
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through a knowledge of g\{m, n). Indeed we have already 
found expressions for g\{m,n) for all members of some 
infinite families of finite groups including all metacyclic 
groups (from which cyclic, dicyclic and dihedral groups 
are specific cases [S^), symmetric groups of permuta¬ 
tions (which include all other finite groups as specific sub¬ 
groups) and alternating groups. Work is in progress to 
find expressions for gx{m,n) for Chevalley groups, which 
are specially interesting because they are the finite coun¬ 
terparts of Lie groups [34] ■ Since all these constructions 
have a purely mathematical character, they shall be pre¬ 
sented elsewhere and we refrain from giving here more 
details. We just give only one example of g\{m,n) that 
describes the infinite family of dihedral groups D 2 q (i.e. 
the group of automorphisms of a regular polygon). These 
groups are noncommutative and have even order X = 2q 
and we find 

g 2 q{m,n) = d, (^0, to-I--I- (23) 

+qd2 ^0, d2 ^0, -b d2 ^0, ^ 


For example, q = 3 describes the group of rotations and 
reflections that fix an equilateral triangle, q = A the ones 
that fix a square, q = 5 a regular pentagon, etc. Setting 
g to a given value, Eq. (231 automatically provides the 
corresponding Cayley table of the group. For example, 
for g = 3 we obtain the Cayley table of Dq, given by the 
operation of ge{m,n) (to is listed in the rows and n in 
the columns) 


ge{m,n) 

0 

1 

2 

3 

4 

5 

0 

0 

1 

2 

3 

4 

5 

1 

1 

2 

0 

4 

5 

3 

2 

2 

0 

1 

5 

3 

4 

CO 

3 

5 

4 

0 

2 

1 

4 

4 

3 

5 

1 

0 

2 

5 

5 

4 

3 

2 

1 

0 


Thus, by using ge^m, n) in Eq. (22), we find the six differ¬ 
ent fractal objects FpQnf{x) obtained from this dihedral 
group. 

Although continuity and differentiability are generally 
lost through the action of the fractal decomposition, the 
fractal objects inherit any symmetry A(a:) of the function 
/(x). Thus if a transformation x —>■ A(x) is such that 
/(A(x)) = f{x) then, we have FpAn/(A(x)) = FpXnf{x) 
as well. For example, if f{x) is periodic, i.e. f{x + L) = 
/(x), so are also the fractal objects necessarily. 

We also note that any permutation of the fractal ob¬ 
jects Fp\nf{x) through the label n, i.e. any bijection that 
sends Fp\nf{x) to Fp\n'f{x), where both n and n' are 
G [0, A — 1] leaves /(x) invariant. Thus, /(x) is invariant 
under the action of the symmetric group of permutations 
of order A, S\, that exchanges the fractal objects in its 
decomposition. 

We have introduced the fractal decomposition in terms 
of the operator Fp^n acting on /. With this in mind, the 


result of the theorem can also be formally expressed as 


A-l 

1 = ^ FpAn (24) 

m—0 

where 1 is the identity operator. In this way no reference 
to a particular / is made. One can easily show that the 
set of A operators Fp^n forms a finite ring under the 
operations 


FpAn * FpAn' — FpA(ra+ra') (25) 

FpAn © FpAn' = FpA(n.n') (26) 

If A is prime this ring becomes a finite Galois field Fa of 
A-characteristic under these operations. 

As in quantum mechanics the operators FpAn can also 
be made to act on vectors. Care must then be taken 
with their proper ordering, since these operators do not 
generally commute. For example, one can define the ex¬ 
ponential operator e = in terms of its 

Taylor series. Trivially, one has then, in general, 

A-l 

^ gFpA„ ^27) 

m—0 


i.e., one cannot factorize such operator because of the 
FpAn being non-commutative in general. This is remi¬ 
niscent of quantum mechanics, where one finds a similar 
situation, leading to the Baker-Campbell-Hausdorff for¬ 
mula and the Trotter-Suzuki decomposition naimssj 
which leads also to fractals m- 

We can then summarize our solution to the problem 
posed in the introduction: It is always possible to bring 
any Hamiltonian to the form of Eq. ([^. One then has 


Hn — FpAnl^ 


(28) 


with FpXnH given by Eq. ( [I7| ). Or, even more generally, 
by Eq. 


IV. EXTENSION TO COMPLEX-VALUED 
FUNCTIONS 


Let z = |z|e*‘'’ be a complex number with g? denoting its 
phase and \z\ its modulus. We have that z = x + iy, with 
X = (z-|-z)/2 and y = {z — z)f{2i) being real numbers (z 
denotes the complex conjugate of z) and it is clear that 
we can extend the definition of the digit function to the 
complex plane as follows 

dp(fc, z) = dp(fc,x)-bzdp(fc,j/) (29) 

The numbers dp(fc,x) and dp{k,y) are non-negative in¬ 
tegers in the interval [0,p — 1] the digit function dp(fc, z) 
is a Gaussian integer in the square with vertices (0, 0), 
{p — 1, 0), (0, p — 1), (p — 1, p — 1) (vertices included). 
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FIG. 4: (Color online) Plot of the analytic complex function f{z) — sin 2 (right) and their associated fractal objects F 2 , 2 , 71 /( 2 ) (i.e. p — 2, X — 2) 
for the values of n = 0 (left) and n — 1 (center). Isochromatic lines denote regions of constant phase. The modulus of f{z) thus varies along 
isochromatic lines. The window displayed is rc G [—5,5] (horizontal axis) and y G [—5,5] (vertical axis). 


Since \z\ > max{|x|, |?/|}, we have 


and we have, as above 


Liogp Id. 


z= (sign(a::)dp(A:, |x|)sign(?/)dp(A:, |y|)) 

(30) 


k — — oo 


We define the complex digit function Cp(k,z) as 

Cp{k, z) = sign{x)dp{k, |x|) + i sign(?/)dp(fc, |y|) (31) 


so that Eq. (301 is most concisely expressed as 
Liogp IdJ 


= Y 


k— — oo 


(32) 


We note that Cp(k,z) yields a Gaussian integer in the 
square with vertices {—p + 1, —p + 1), [p — 1, —p + 1), 
(—p + 1, p — 1), and (p — 1, p — 1) (vertices included). 
Since we also have 


z = z = e 


Liogp IdJ 

£ /dp(fc,|z|) 

k— — oo 


(33) 


p-i 

f{z) = Y FpA„/(^) 


n—0 


(36) 


and the fractal objects Fp^^ form again a Galois field 
if A is prime. As above, the factor d;),(0, |fc| + n) can be 
replaced in Eq. (35) by an operator gx(k,n) with the 
latin square property and Eq. (36) still holds (although 
the Galois field property is then generally lost and we 
have then a finite ring). 

In Fig. (|^ the complex-valued fractal objects 
F 2 , 2 ,nSinz (i.e. p = 2, X = 2 and f(z) = sinz) cal¬ 
culated from Eq. (35) for the values of n = 0 (left) and 
n = 1 (center) are shown. The sum of these objects 
yields f{z) = sin 2 (right), as calculated from Eq. (36). 
We emphasize the fact that although the fractal objects 
generally fail to be analytic, the symmetries of f{z) are 
also symmetries of the fractal objects. Thus, if for exam¬ 
ple f{z) is invariant under a Mobius transformation so 
that 


the following identity holds 

Cp{k, z) = e*‘^dp(fc, Izl) = T^dp(fc, \z\) (34) 

With all these definitions and straightforward identi¬ 
ties in mind, the fractal decomposition can now be easily 
extended to complex numbers in an analogous manner 
to the one described above for the reals. The pXn fractal 
decomposition of a complex function f{z), f : C —)■ C, 
takes the form 

FpXnfiz) = (35) 

^ Liogpl/G)IJ 

Y P''cp{kj{z))dx{0,\k\+n) 

k— — oo 


for some integers a, b, c, d, such that ad — be > 0, then 
each fractal object in the decomposition is also invariant 
under the same Mobius transformation. In the case of 
f{z) = sinz in the figure, since the function exhibits a 
periodic behavior, the fractal objects have this overall pe¬ 
riodicity as well and are somehow ‘phase locked’ to f{z) 
(albeit with a more detailed and fine structure). This is 
clearly noticed in the figure: Stripes with a certain color 
(the isochromatic regions denote constant phase regions) 
occupy the same positions in the fractal objects as in 
f{z). Also, the zeros of f{z) are zeros of the fractal ob¬ 
jects as well (although the contrary is not true). 


A(A-l) 











V. THE GENERALIZED p\n FRACTAL 
DECOMPOSITION 


A wide family of fractal decompositions, that include 
the FpAn decomposition as a specihc case, can now be 
constructed by further inspecting the structure of Eq. 
(15) which is instrumental to establish the decomposi¬ 
tion, Eq. ( [T8| ). We state this result rigorously as a the¬ 
orem, prove it, and give further concrete examples of 
fractal decompositions. 

Theorem. Let S he the set of integers G [0, A — 1] 
(A > 2, A S Nj and let the digit function d\{k,x) for 
a; € K and k G Z be defined by Eg. Let also 

g\ (m, n) : Z X Z ^ S be a surjective application such 
that gx{m,n) = to), d>,(0, n)) is also a law of 

composition (magma) S x S ^ S with the latin square 
property. Then, for any injective function h : S ^ M"*" 
(where denote the non-negative real numbers) the fol¬ 
lowing relationship holds 


A-l 


A-l 


'^h{gx{m,n)) = '^h{j) = ax 


(38) 


integers G S. Each different to leads to a different per¬ 
mutation h(j') of the same distinct A values. Since all 
these permutations contain the same distinct elements in 
the image set, the sum of the latter is the same, inde¬ 
pendently of TO. This proves Eq. (38). From the latter. 


since all h{j) are distinct and non-negative, cta ^ 0 and 
we have 


— ^h{gxi\k\,n)) = 1 


A-l 


(41) 


for any k G Z. We now observe that, since each 


h {gx (I/cI, n)) is non-negative. 


< 1 and. 


n—0 


j=0 


hence, for any well dehned function f{x) expressed in 

radix p as f{x) = sign/(a:) l/(a;)l) 

the fractal object GpXnf{x) given by Eq. (39) is neces¬ 
sarily well defined as well, being Cauchy-convergent on 
each point x to a certain real value (even when continuity 
and differentiability, which both depend on the specihc 
forms of the functions h and /(x), is not generally war¬ 
ranted). Eurthermore, 


where ax G K’*' is a nonzero constant, independent ofm. 
Furthermore, for any function f{x) : K —> M, there exist 
A distinct fractal objects GpXn defined as 

GpXnfix) = (39) 

, . Llogp l/(a:)IJ 

^sig^ ^ pMp(fc,|/(x)|)h( 5 A(|fc|,n)) 

that constitute a generalized pXn fractal decomposition, 
satisfying the following relationship 

A-l 

/(x) = ^ GpXnfix) (40) 

n=0 


Proof: The surjective application gxim,n) of any two 
integers m G Z and n G Z homomorphically maps 
the cartesian product Z x Z to the cartesian product 
S X S, where S is the set of integers in the interval 
[0, A — 1] (A G N), and returns an integer also in S, since 
gx{m,n) = 5>(d;),(0 ,to), d>(0,n)) G S. Now, since the 
law of composition (;A(dA(0, to), dA(0, n)) has the latin 
square property, all A distinct integers j G S, 0 < j < 
A — 1 appear in each row dA(0,TO) of its Cayley table. 
The injective function h : S —>■ K’*' maps each of these 
distinct A integers to a distinct non-negative real num¬ 
ber h (pA(dA(0, to), dA(0, n))) (distinctness is preserved 
because the function h is injective). Therefore, the latin 
square (/A(dA(0, to), dA(0, n)) with entries in S becomes a 
latin square h (pA(dA(0, to), dA(0, n))) with entries in K+. 
If we now hx to and change n so that dA(0,n) changes 
from 0 to A — 1 it is clear that h ((/A(dA(0, to), dA(0, n))) 
ranges over all A distinct non-negative images /i(j) of all 
distinct A entries j G S which are a permutation of all 


A-l A-l . J., s 

r: : Sign/ x) 

2^ GpXnfix) = 2_^ - X 


n=0 


CTa 


(42) 


n=0 

[logp l/(a:)IJ 

X P^<ipik,\fix)\)higxi\k\,n)) 

k——oo 

[logp l/(®)IJ 

= sign/(x) p^dpik,\fix)\) X 


k— — c 


\ X ra=0 / 


[logp l/(a;)|J 

= sign/(x) Y p''dp(A:, |/(x)|) =/(x) 

k— — oo 


where Eq. (41) has been used. This completes the proof 
of the theorem. □ 


The theorem also holds for a function of several real or 
complex variables /(z, w,...) because no property of fix) 
is used in proving the theorem other than the fact that it 
admits a radix-p expansion. Eunctions with an arbitrary 
number of variables also generally admit such expansion. 
Also, any complex-valued function fiz,w,...) of an ar¬ 
bitrary number of variables, admits a radix-p expansion 
by means of the complex digit function Cp(fc, /(z, w ,...)). 
In this latter case one has 


GpXnfiz,W, ...) = 


(43) 


1 

cr\ 


[logp 

Y P''cpik,fiz,w,...))higxi\k\,n)) 


k— — oo 
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We then have, 

A-l 

f{z, w,...) = Y^ GpXnfiz, w,...) 

n—0 

Corollary. Let Ao,Ai,...,Am be natural numbers, all 
> 2 and let f be an arbitrary function that admits a 
radix-p expansion. Then, 

Ai —IA 2 —1 Am ~1 

f— A ' A ' A ' GpAjMMM ■ • ■ GpA2n2 GrpAini / 

no=On2=0 TiM—O 

(44) 

Proof: Since / admits a radix-p expansion, so does 
it fi = GpAim/ as well, as a result of the theo¬ 
rem, and so does f 2 = GpX 2 n 2 fi — GpX 2 n 2 GpXinif' 
One then has that fM = GpA„n„/M-i = 
Gpx^timGpXm -1 ''' GpXini f is also well defined. 
We now note 


Am —1 Am —1 

^ ^ /m = ^ ^ Gt^Am^^m/m —1 = /m —1 

riM—O riM—O 


(45) 


from Eq. (40). Therefore 


Am-1 —1 


Am —1 — 1 Am — 1 

AZ ^ = AZ fM-l = fM-2 

nM-l=0nAf=0 nM-l=0 


(46) 


also from Eq. (40) and we have, as well (by induction) 

Ai — 1 Am —1 1 Am — 1 Aq — 1 

E--- E E /m=E/i = / ( 47 ) 

ni=0 nM-i=0nM=0 ni=0 


which is Eq. (441. □ 


We now give some examples of the theorem. In all 
cases h(j) constitutes an example of a suitable injective 
function. 

Example 1: If h{j) = j (j € S'), and qx(rn,n) = 
dA(0,m-|-n) we have cjx = G^xn 

given by Eq. (391 coincides with FpA™ giveiiby Eq. ([T7[). 


Example 2: If h{j) = In(I-l-j), with j G S and gx(jn, n) = 
dA(0,m -I- n) we have <Jx = In A! in Eq. (38) and GpAn 
given by Eq. (39) takes the form 


r. _ signfix) ^ 

GrpAn/((^) — , . . X 


(48) 


In A! 

Llogp l/(a:)IJ 

X AA p''dp(fc,|/(x)|)ln(l + dA(0, |fc|-hn)) 


= —OO 


Example 3: If h{j) = j^, with j G S and q G N and 
gx{m,n) = dA(0,m -I- n) we have, by the Seki-Bernoulli 
formula that 


Sw 9 + 1-^ ' 


(49) 


in Eq. (38). Here are the Bernoulli numbers Bq = 1, 
Bi = i, B 2 = i, B 3 = 0, Bi = -i, etc. Then GpA„ 
given by Eq. ([3^, takes the form 


GpXnf (^) — 


sign/(a;) 


(q\ (A-i)‘i+i-^ 


(50) 


e;.o(?) 

Llogp |/(a;)|J 

X E /dp(A:, |/(a:)|) [dA(0, |A:|-hn)]”^ 

k— — oo 

Let us apply this to the case A = 4 and 9 = 3 to better 
illustrate this example and the proof of the theorem. We 


latin square 54(771 

77) 

= 

d 4(0 

d 4 ( 0 , 771 - 1 - 77 ) 

0 

1 

2 

3 

0 

0 

1 

2 

3 

1 

1 

2 

3 

0 

2 

2 

3 

0 

1 

3 

3 

0 

1 

2 


with m listed in the rows and n in the columns, for values 
in the set S. Now, this latin square is sent by the injective 
application h = to another latin square 


[d4(0,777 -1- 77 )]^ 

0 

1 

2 

3 

0 

0 

1 

8 

27 

1 

1 

8 

27 

0 

2 

8 

27 

0 

1 

3 

27 

0 

1 

8 


The sum of any row or any column of this latin square is 
a 4 = 36 as predicted by the Seki-Bernoulli formula, Eq. 
(|4^, for A = 4 and q = 3, 


(74 = E 

e=o 

^3\^ 
vOy 4 
= 36 


3\ 3 


4-£ 


4- 




3\ 33 1 

ijy 2 


3A(^1 
2/2 6 




Therefore, Eq. ( [^ reduces in this specific case to 
5 _ sign/(a:) 

Gp4nf\X) — 36 ^ (® 7 ) 

Llogp |/(a;)|J 

X E /dp(fc, |/(x)|)[d 4 ( 0 , |fc|-hn)]^ 

= —00 


Example 4-' If h{j) = r-l, with j G S and r > 0, r S K"*" 
and gx{m,n) = dA(0,m -I- n) we have a geometric series 
and, hence, 


crx = 


pA _ 1 

r — 1 


(52) 
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in Eq. (38). Therefore, from Eq. (39), 


^pXnf (^) — 


_ (r - l)sign/(x) 


— 1 


(53) 


[logp l/( 2 :)IJ 


k——oo 


change in the statistical system. Indeed, from the asymp¬ 
totic limit p —>■ oo of Eq. (54), (cf. Eqs. ([^ and (|^) 


/ = lim dp{-k,f)p 

p—^oo 

fc=-[logp /J 

= doo(0,/) 


-k 


(62) 


VI. CONNECTION TO STATISTICAL 
MECHANICS 

A formal connection can be established between digit 
expansions in positional number systems and the statisti¬ 
cal mechanics of systems with a discrete energy spectrum. 
In this section we sketch this connection. We focus on 
non-negative real valued functions f{x,y,...) > 0. The 
radix-p expansion of such a function is given by Eq. (§ 
as 


/ 


[logp /J 


k— — oo 


Y p ^dp(-fc,/) 

fc=-Liog /J 


(54) 


Since / > 0 we can formally define a ‘Helmholtz free en¬ 
ergy’ F such that / = with /3 = 1 /(A:bT) denoting 

a Lagrange parameter, with T the ‘temperature’ and ks 
being ‘Boltzmann’s constant’. We can then consider a 
‘statistical system’ with a discrete ‘energy’ spectrum de¬ 
scribed by a ‘canonical ensemble’ and where E denotes 
energy, Eq the energy of the ground state, e a ‘character¬ 
istic energy’ (governing the spacing of the energy levels), 
and fl{E, n, v) the number of microstates with energy E, 
number of particles n, volume v (other extensive variables 
that may also be present). If we introduce the following 
correspondence 


n 

= X V = y etc. 

(55) 


= fix,y) 

(56) 

E 

e 

= k 

(57) 

Eq 

e 

= -[^ogpfix,y)\ 

(58) 

e 

keT 

= /3e = hip 

(59) 

,n,v) 

= <ip{-k,f{x,y)) 

(60) 


then Eq. ([54| becomes 


we find that / equals its own most signihcative digit 
when written in an infinitely large radix (i.e. there is 
only one digit in this representation: the most represen¬ 
tative one). This limit, because of Eq. (59), is similar 
to taking in the ‘statistical system’ the limit T —>■ 0, so 
that the system is fully described by the ground energy 
state Eq. a simple example may help to clarify this. 
The positive number / = = 204 (e.g.) has only 

one digit in radix p > 205 (in radix p = 205 one has 
204 = 204 • 205° and hence, only one power of 205 is 
needed to describe the number). In the decimal radix 
p = 10, we have that 204 = 2 • 10^ -I- 0 • 10^ -I- 4 • 10° 
has three digits, 2, 0 and 4, and three powers of ten 
are needed. Finally, 204 in radix p = 2 has 8 digits, since 
204 = 1-2'^+ 1-2^+ 0-2^+ 0-2'^+ 1-2^+ 1-2'^+ 0-2^+0-2^. 
Therefore, if we lower the radix we need more digits to de¬ 
scribe the same number. The statistical correspondence 
established above shows, through Eq. (59), that lowering 
the radix p is equivalent to increasing the temperature T 
in the statistical system at constant Z = e~^^: In the 
low temperature limit, the system is found in only one 
energy level, i.e. the ground state; as temperature is in¬ 
creased, excited energy levels become occupied and more 
energy levels are needed to describe the state of the sys¬ 
tem. This increased number of energy levels of the statis¬ 
tical system mirrors the increased number of digits needed 
at lower radix to represent the number . 

We should now clarify how the pXn fractal decompo¬ 
sition enters in this discussion. To understand this let 
us consider the possibility that the number of particles 
in the system n can fluctuate so that n can range from 
0 to a maximum number N^ax of particles. The sys¬ 
tem is in contact with a reservoir at chemical potential 
/i with which it exchanges particles. A suitable, natural 
description of this system is provided in statistical ther¬ 
modynamics by the grand-canonical potential S . We 
have 


e 


-/3H 



n—0 


^-PF{n)^l3im 


(63) 


e = 'Y ^{E,n,v)e (61) 

E—Eq 

which indeed constitutes the definition of the canonical 
partition function Z of the statistical system. It is then 
interesting to note that the change of radix p in the num¬ 
ber theoretical description of a non-negative real num¬ 
ber mathematically mirrors the effect of a temperature 


where F{n) is the free energy of a cluster of n parti¬ 
cles. This expression defines the grand-canonical parti¬ 
tion function. From here we observe that if we identify 
Nmax with A — 1 where A is the number entering in a 
generalized pXn decomposition, we then note that we can 
further identify 

e-PF{n)^Ppu ^ GpA„e-^= 


(64) 
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and then Eq. (63) becomes 


A-l 


= Eg 


pXn^ 


(65) 


n—0 


which is consistent with Eq. (40). We then note that the 


probability V{n) of observing a cluster with n particles 
is given by 


V{n) = 


g-0F(n)^j3nn 


= 


pXn^ 


( 66 ) 


This expression is remarkable, since it is valid also if 
is complex-valued. It yields, in any case, a real-valued 
probability. To see this note that if we make the transfor¬ 
mation /3S —>■ /3S -|- ia (with a an arbitrary real number) 
we have, because of Eqs. (|M| and Eq. (43), 


so that, again, from Eq. (69) the consistent relationship 

^A-l 

=o' 


e = X)Tt=o GpA'A/e is warranted. The generalized 


ensemble describes small systems out of the thermody¬ 
namic limit [55MD] and has been successfully used by 
Chamberlin to derive a model for ferromagnetism |39j 
and a mean-field mesoscopic Landau theory [40]. We 
note that the subdivision potential £ is equal to zero in 
the thermodynamic limit as a consequence of the Gibbs- 
Duhem relationship, which precludes that all extensive 
variables in the system simultaneously fluctuate. Follow¬ 
ing previous work by Vives and Planes m we found in 
that the subdivision potential £ of Hill’s Nanother¬ 
modynamics can be connected with the Tsallis entropic 
form [5] 


Oq — 

l-q 


(71) 


^pXn 


—p£^—ia 


) = e-*“GpA„e- 


(67) 


Therefore, by introducing the same transformation ^ —>■ 
S-|-ia in Eq. (66) we see that the latter remains invariant 


V{n) = e'5=+*“GpAne-'^' 


— p£l — lOL 


= 


pXn^ 


( 68 ) 


Here g G K is the entropic parameter and pj is the prob¬ 
ability of a microstate j. The interest of Tsallis entropy 
lies in its connection with (multi) fractals and complex 
systems [5|imi33| being a most prominent form of super¬ 
statistics [23|. The connection found in [32] establishes 
that 


We have thus proved the phase invariance of Eq. (66) 


and, thus, our methods are also able to describe Hilbert 


spaces, with Eq. (66) giving the observable probability 


(Born rule) of a quantum state governed by a ‘wave func¬ 
tion’ where /IS can be an arbitrary complex num¬ 

ber. Indeed, in quantum mechanics /IS can be thought 
as having complex part equal to iS/h where S is the La- 
grangian action and h the Planck constant. We observe 
also the remarkable fact that the low-temperature limit 
/I —>■ oo of a statistical system admits a dual interpre¬ 
tation under this connection as the limit 5'//i —>■ oo of a 
quantum system. Thus, if we take the radix p = \_S/h\ = 
ry —>■ 00 m the ’ground state’ of the statistical system is 
nothing but the classical limit of a physical system where 
no nontrivial partitions are observed (note that in this 
regime, Eq. ^ holds). In fact, we explained in [57] how, 
when demanding that the radix rj has the least economy, 
physical laws can be derived in a unified manner. 

Coming back to our connection to statistical mechan¬ 
ics, If we now further consider a situation where the vol¬ 
ume of the system V = Afv can also fluctuate in multiples 
Af of an specific volume v at the pressure of the reservoir 
p we now have that this situation is well defined by a 
generalized free energy (subdivision potential) £ 
as 


e-f}£ ^ g g-/3H(AA)g-/3p«N (gg) 

J\f=0 

This relationship defines the partition function for the 
generalized ensemble. We find, again, by introducing 
A' = Afmax that we can identify 

= GpA'A^e-^^ (70) 


The entropic parameter q controls the ‘distance’ to the 
thermodynamic limit, which is approached when q —> 1. 
In that limit, by using L’Hopital’s rule we observe from 
this expression that, indeed, £ = 0, as mentioned above. 


VII. CONCLUSIONS 

In this article we have established a general mathemat¬ 
ical mechanism to partition any quantity into a finite set 
of fractal functions whose ordinary sum is equal every¬ 
where to the original quantity. We have shown that there 
are infinite families of possible partitions indeed, all them 
exact, and we have explained how these partitions can be 
systematically constructed. These results have also been 
extended to complex-valued functions. 

The decomposition thus allows a finite family of fractal 
objects to be derived from an arbitrary ‘mother’ function 
f(z), which can be continuous and analytic. Perhaps a 
most interesting aspect of the construction is that, al¬ 
though the fractal objects generally fail to be continuous 
and analytic, they possess all internal symmetries of f{z). 

A formal connection between our mathematical ap¬ 
proach and statistical mechanics has also been briefly 
sketched. The interest of this connection, which shall 
be further explored elsewhere in the context of specific 
applications, lies in the fact that it provides a number- 
theoretic foundation for statistical mechanics. We have 
shown that digit expansions in positional number systems 
are naturally connected to partition functions as the ones 
arising in statistical mechanics. Interestingly, the radix 
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change in the representation of numbers mirrors the ef¬ 
fect of a temperature change in a statistical system. The 
connection also constitutes a pathway to the investiga¬ 
tion of some certain contemporary approaches to statis¬ 
tical thermodynamics in which fractality, complexity and 
cooperativity play a central role 0 138114(71143ll4^ . 

The radix representation of numbers is a topic that 
seems to have been neglected in Physics because of the 
obvious fact that a number is ‘always the same’ indepen¬ 
dently of the radix used. Nonetheless, we have reasoned 
here why we do believe that the radix of numbers has an 
important physical meaning. Although a number is itself 
invariant under a radix change, its specific partitions in 
powers of the radix are not: A number, as represented in 
a standard positional number system, is also a partition 
function, where digits accompany powers of the radix, the 
total number of such powers being dependent both on the 
number and on the radix. Therefore, the change of the 
radix of physical numbers has a physical impact because 
it governs the structure of the partitions of the number 
and the latter have a major importance in physics when 
one considers the quantum regime, stochastic phenom¬ 
ena, or the behavior of systems with many degrees of 
freedom. In a recent article HZ], we have identihed a 
physical radix that governs the transition between clas¬ 


sical and quantum mechanics and which is related to the 
Lagrangian action. When this radix is optimal, in a sense 
that it minimizes the size of a certain matrix called radix 
economy (or digit capacity), physical laws are derived in 
a straightforward manner and the principle of least ac¬ 
tion of classical mechanics is encompassed as well as the 
Hilbert space of quantum mechanics. In m we have 
also established a relationship between the optimal radix 
(given by the Lagrangian action) and Boltzmann entropy 
(see also [50]) and in the present paper we have sketched 
how a change of the radix is related to a change of tem¬ 
perature in a canonical ensemble. 

We believe that the approach presented here makes 
fractals amenable to theoretical investigation from a new 
perspective, in connection with problems involving quan¬ 
tization [TM1I35|. In quantum mechanics, symmetries 
are related to conservation laws that dictate which kind 
of physical processes at a fundamental level are possi¬ 
ble, constraining them. Such symmetries and constraints 
are subtly and beautifully related through group the¬ 
ory. This article also establishes a general framework 
that connects fractals and finite group theory in a way 
that seems natural, and we hope that it may prove a use¬ 
ful tool for applications in both statistical mechanics and 
quantum field theory. 
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